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Abstract 

Background: Upper tract urothelial carcinomas (UT-UC) can invade the pelvicalyceal system making differential diagnosis of 
the various histologically distinct renal cell carcinoma (RCC) subtypes and UT-UC, difficult. Correct diagnosis is critical for 
determining appropriate surgery and post-surgical treatments. We aimed to identify microRNA (miRNA) signatures that can 
accurately distinguish the most prevalent RCC subtypes and UT-UC form the normal kidney. 

Methods and Findings: miRNA profiling was performed on FFPE tissue sections from RCC and UT-UC and normal kidney 
and 434 miRNAs were significantly deregulated in cancerous vs. the normal tissue. Hierarchical clustering distinguished UT- 
UCs from RCCs and classified the various RCC subtypes among them. qRT-PCR validated the deregulated expression profile 
for the majority of the miRNAs and ROC analysis revealed their capability to discriminate between tumour and normal 
kidney. An independent cohort of freshly frozen RCC and UT-UC samples was used to validate the deregulated miRNAs with 
the best discriminatory ability {AUC>0.8, p<0.001). IVlany of them were located within cytogenetic regions that were 
previously reported to be significantly aberrated. miRNA targets were predicted using the miRWalk algorithm and ingenuity 
pathway analysis identified the canonical pathways and curated networks of the deregulated miRNAs. Using the miRWalk 
algorithm, we further identified the top anti-correlated mRNA/miRNA pairs, between the deregulated miRNAs from our 
study and the top co-deregulated mRNAs among 5 independent ccRCC GEO datasets. The AB8/13 undifferentiated 
podocyte cells were used for functional assays using luciferase reporter constructs and the developmental transcription 
factor TFCP2L1 was proved to be a true target of miR-489, which was the second most upregulated miRNA in ccRCC. 

Conclusions: \Ne identified novel miRNAs specific for each RCC subtype and UT-UC, we investigated their putative targets, 
the networks and pathways in which they participate and we functionally verified the true targets of the top deregulated 
miRNAs. 
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Introduction 

Renal cell carcinoma (RCC) represents 2—3% of all cancers and 
accounts for approximately 90% of all kidney malignancies. Apart 
from surgery, it is both chemotherapy and radiotherapy resistant 
and it is composed of various morphologically and cytogeneticaUy 
distinct subtypes. The most prevalent subtypes are clear cell RCC 
(ccRCC, 75-80%), papiEary RCC (papRCC, 10-15%) and 
chromophobe RCC (chRCC, 5%) [1]. Distinguishing RCC 
subtypes is of clinical importance because they have difiFerent 
prognoses and consequendy different management schemes [2]. 
However, morphology-based distinction is not always conclusive 
since some subtypes may have overlapping or related morphologic 



features. For the purpose of targeted therapy it is especially 
important to classify the different subtypes of RCC. The 
histological types arise from different cells of origin in the kidney, 
different constellations of genetic alterations [3], and expression or 
mutation in different oncogenic pathways. Therefore, different 
subtypes offer different molecular candidates for targeted therapy, 
such as Tyrosine Kinase Inhibitors, Sorafenib and Sunitinib, 
mTOR inhibitors, Everolimus and Temsirolinius, etc. There is 
growing evidence that variability in response rates may be linked 
to sub-classification [4] . Therefore, new biomarkers are needed in 
order to improve the identification and diagnosis of renal tumor 
subtypes. Recent data suggest that RCC classification through 
microRNA (miRNA) expression profiles is highly accurate [5,6]. 
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Upper tract urothelial carcinoma (UT-UC) is a relatively 
uncommon form of kidney cancer arising from the urothelial 
lining of the renal pelvis and calyces. UT-UC accounts for the 
majority of bladder cancer; however, it only accounts for about 
7 % of renal neoplasms [7] . UC of the renal pelvis is an aggressive 
tumour, which may invade the renal parenchyma, mimicking 
primary renal cell carcinoma. Similarly, advanced RCC can 
invade the pelvicalyceal system. This can make differential 
diagnosis of RCC and urothelial carcinoma of the renal pelvis 
difScult. Correct diagnosis is critical for determining appropriate 
surgery and post-surgical treatments. For instance, UT-UC 
including renal pelvis, calyces and ureters will require radical 
nephrectomy with ureterectomy and bladder cuff resection. 
However, RCC will require only partial or radical nephrectomy 
without extensive ureter resection. Correct diagnosis is critical for 
determining appropriate surgery and post-surgical treatments. 
Therefore, it is of major importance to identify biomarkers that 
can accurately distinguish UT-UC from RCC [8,9]. 

MicroRNAs (miRNAs) are small non-coding RNAs of approx- 
imately 19-23 nt size, shown to regulate gene expression at the 
post-transcriptional level, by binding through partial sequence 
homology to the 3' UTR of mammalian target mRNAs and 
causing translational inhibition and/or mRNA degradation [10]. 
miRNAs have been identified as key regulators in many biological 
processes including development, differentiation, apoptosis and 
proliferation. They are aberrandy expressed or mutated in 
cancers, suggesting that they may play a role as a novel class of 
oncogenes or tumour suppressor genes [1 1]. Deregulated miRNAs 
have been found in many tumours, including kidney cancer and 
carcinomas of the urinary tract [5,12,13,14,15,16,17,18,19,20,21]; 
however, not much is known regarding their expression profile in 
upper tract urothelial carcinomas. Furthermore, in RCC the 
results are not always consistent [12,13,17,18,22,23,24,25]. 
Additionally, many deregulated miRNAs are located in chromo- 
somal regions that are frequently altered in cancerous situations 
[26]. 

To our knowledge there are no accurate biomarkers for early 
detection of both RCC and UT-UC, thereby accounting for late 
diagnosis and subsequent poor prognosis. Better understanding of 
RCC and UT-UC pathogenesis is required for identifying 
biomarkers and developing novel targeted therapies. In the 
present study, we examined the miRNA expression profdes 
mainly of ccRCC, as well as in a smaller cohort of papRCC, 
chRCC and UT-UC samples. We identified novel miRNAs that 
are specific for each tumour entity, we investigated their putative 
targets, and we explored the networks and pathways in which 
these miRNAs participate. Finally, we functionally verified the 
target gene of the top deregulated miRNA in ccRCC. 

Materials and Methods 

Patient samples 

Formalin-fixed and paraffin embedded (FFPE) tissue samples 
from patients with kidney cancer and upper tract urothelial 
carcinoma were retrieved from the archives of the Asklipieio 
General Hospital, Athens, following research ethics board 
approval from the same Hospital. Paraffin sections from each 
specimen were reviewed by two independent pathologists and 
histologically classified according to the WHO classification and 
staged according to the UICC-TNM classification (2009). Nuclear 
grade was scored according to the Fuhrman classification system. 
Kidney specimens were collected between 2007 and 201 1. Written 
informed consent was obtained from all patients included in this 
study. The study protocol was approved by the Ethics Committees 



of the Asklipieio General Hospital, Athens and the University of 
Cyprus. Twenty-four ccRCC, 5 papRCC, 3 chRCC, 5 UT-UC 
patients and 40 normal kidney tissue (parenchyma/ cortex) samples 
which were used as controls, were enrolled in the present study 
(Table 1). None of the patients had a history of previous bladder 
cancer or concomitant bladder cancer at UT-UC diagnosis. The 
criteria to select normal kidney tissues were that the normal kidney 
cortex tissue was from the same individual, from the same kidney, 
from the same surgical specimen, and collected under the same 
standard operation procedure conditions: if possible the sampling 
of normal cortex was done leaving ~2 cm of macroscopicaUy 
healthy cortex adjacent to the tumor untouched. ccRCC are 
macroscopicaUy characterized by a sharp demarcation (fibrous 
pseudocapsule) between tumor and surrounding non-tumoral 
tissue. This is also reflected microscopically, as RCC are sharply 
demarcated from the surrounding tissue. The non-tumoral tissue 
was therefore obtained in macroscopicaUy clear-cut benign tissue. 
Haematoxylin and eosin (H&E) staining and immunohistochem- 
istry (IHC) for several markers (AE1/AE3 keratins, CK7, CK8, 
CK18, CK19, EMA, CDIO, Vimentin, EGFR, Ki67, NSE, p53 
and S-lOO) was performed for all RCC samples. Anti-vimentin 
staining was also performed for UT-UC samples (Figure 1). 

RNA extraction and microRNA microarray analysis 

Total RNA containing small RNAs was isolated from 18 
ccRCC, 5 papRCC, 3 chRCC, 4 UT-UC and 20 normal tissue 
samples using the High Pure FFPE RNA Micro Kit (Roche 
Applied Science, IN). The quality of the total RNA was verified by 
an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, 
CA) and the concentration was measured spectrophotometrically 
(Nanodrop technologies, Montchanin, DE). Tissue sections were 
deparaffinized with xylene and ethanol washes, treated with 
protease and total RNA containing small RNAs was isolated using 
the High Pure FFPE RNA Micro Kit (Roche Applied Science). 
Total RNA (0.5 |ig) from each sample and reference was labeled 
with Hy3 fluorescent label, using the miRCURY LNA microRNA 
Hi- Power Labeling Kit (Exiqon, Vedbaek, Denmark), according to 
the manufacturer protocol. The Hy3-labeled samples were 
hybridized to the miRCURY LNA microRNA Array, 6* gen 
(Exiqon, Vedbaek, Denmark), containing capture probes that 
target aU miRNAs for all species registered in the miRBase vl8.0. 
The hybridization was performed using an AgUent hybridization 
SureHyb chamber and gasket slide kits. After hybridization, the 
microarray slides were scanned at 10 |im using the High- 
Resolution Microarray Scanner (Agilent Technologies) and stored 
in an ozone free environment. The image analysis was carried out 
using the ImaGene 8.0 software (BioDiscovery, Inc., USA). 
FUtering was performed based on the signal intensity. Spots with 
no signal above the background (flags 1 and 2) were detected and 
removed. Background correction was performed by subtracting 
the median global background from the median local background 
from the signal intensity. Normalization was performed using the 
quantile normalization algorithm. Normalized data were further 
extracted, pre-processed and sorted with Microsoft Excel®. For 
further data analysis the Matlab® (The Mathworks Inc.) comput- 
ing emdronment was used. MicroRNAs were considered to be 
sigiiificaiitiy differentially expressed (DE) if they obtained a p- 
value<0.05 and a FDRSO.05. Two-way average-linkage hierar- 
chical clustering (HCl) with Euclidian distance was performed with 
Matiab® and Genesis 1.7.2 (Technische Universitaet-Graz, 
Austria) software. Our MIAME compliant microarray data were 
deposited in the Gene Expression Omnibus (GEO; www.ncbi.nkn. 
nih.gov/geo) under accession number GSE48008. 
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Table 1. Cllnicopathological information of the patients. 







ccRCC (n = 24) 


papRCC (n = 5) 


chRCC (n = 3) 


UT-UC (n = 5) 


Median age (range), y 


70 (43-79) 


67 (53-80) 


68 (60-77) 


76 (76-78) 


Gender (male/female) 


20/4 


4/1 


0/3 


4/1 


Stage 


Tl 


13 


4 


0 


3 


T2 


6 


1 


3 


1 


T3 


5 


0 


0 


1 


T4 


0 


0 


0 


0 


Metastatic tissue samples 


2 


0 


0 


0 


Grade 


1 


0 


0 


0 


0 


2 


5 (Tla), 3 (Tib), 2 (T2a), 1 (T2), 1 
(T3a) 


2 (Tla) 


1 (T2b) 


1 (Tl+in situ) 


3 


3 (Tla), 2 (Tib), 1 (T2a), 3 (T3a) 


2 (Tib), 1 (T2a) 


2 (T2a) 


1 (Tl), 1 (T2), 1 (T3) 


4 


1 (T2a), 1 (T2b), 1 (T3a) 


0 


0 


0 


median tumour size (range), cm 


Grade 1 . _ _ _ 


Grade 2 


4 (2-8.5) 


2.75 (2.5-3) 


1 3.2 (0) 


1.5 (0) 


Grade 3 


6.5 (1.5-9) 


5.5 (4.5-8.5) 


9.25 (8.5-10) 


0.7 (0.5-1.5) 


Grade 4 


8.4 (6.5-13) 








Nodal status. 


positive/negative/unknown 


0/24/0 


0/5/0 


0/3/0 


0/5/0 


Metastatic status. 


positive/negative 


2/22 


0/5 


0/3 


0/5 



doi:1 0.1 371 /journal.pone.0091 646.t001 



Quantitative real time polymerase chain reaction (qRT- 
PCR) validation in FFPE samples and an independent 
cohort of freshly frozen tissue samples 

qRT-PGR for the top DE miRNAs was performed for aU FFPE 
samples for which total RNA was initially extracted (24 ccRCCs, 5 



papRCCs, 3 chRCCs, 5 UT-UCs and 40 normal kidney tissue 
samples) including those used for the microarray analysis. 

Additionally, a blinded independent validation set of 40 freshly 
frozen kidney tissue samples (10 ccRCCs, 3 papRCCs, 4 chRCCs, 
3 UT-UCs and 20 normal kidney cortex samples) was used to 
study the assay performance. The validation set included samples 



ccRCC,VIM+ papRCCVIM^ chRCC.VIM- normal MOney, VIM- UT-UC, VIM- 



H&E 



IHC 









■ I ' ■ 















Figure 1. Haematoxylin and eosin (H&E) staining and immunohistochemistry (IHC). Upper pannel: Representative H&E staining from 
ccRCC, papRCC, chRCC and normal kidney tissue. About 5-10 serial tissue sections of 10 |^m were cut from each FFPE block, deparaffinized with 
xylene, hydrated and stained with H&E before microscopic examination. When the proportion of tumour cells was >70% the FFPE block was 
subjected to total RNA extraction. Lower pannel: IHC of FFPE tissue sections using antl-vimentin as primary antibody. VImentIn was predominantly 
seen In ccRCC and papRCC (—70% and —50%, respectively), but only rarely In chRCC (4%) and absent In the normal kidney. VImentIn was also down- 
regulated In the majority of UT-UC cases. 
doi:1 0.1 371/Journal.pone.0091 646.g001 



PLOS ONE I www.plosone.org 



3 



March 2014 | Volume 9 | Issue 3 | e91646 



miRNA Classification of Kidney Tumours 




Figure 2. microRNA profiling. Four-hundred and thirty-four miRNAs were statistically significantly deregulated in all RCC subtypes and UT-UC 
versus the normal kidney. (A) Q-Q (quantile-quantile) plot. Red circles indicate the significantly deregulated miRNAs. (B) Frequencies of the t-scores 
and p-values. The deregulated miRNAs had a p<0.05. (C) The volcano-plot depicts the 434 statistically significantly deregulated miRNAs in ccRCC, 
papRCC, chRCC and UT-UC versus the normal kidney, of which the majority was significantly down-regulated in the cancerous tissue compared to the 
latter. (D) FDR diagram depicting the percentage of FDR with respect to p-value along with a plot of the estimated a priori probability that the null 
hypothesis Tt(0), is true versus the tuning parameter, lambda, X, with a cubic polynomial fitting curve. 
doi:10.1371/journal.pone.0091646.g002 



that met the same .selection criteria with the initial FFPE coliort, 
such as >70% tumor cells in each sample and concordance in 
diagnosis between two independent pathologists. Samples were 
collected at the Asklipieio General Hospital, Athens. Only samples 
for which both pathologists agreed in the diagnosis were included 
in the validation cohort. 

Total RNA was extracted by homogenizing approximately 
0.2 g of frozen tissue with a tissue homogenizer, followed by 
isolation using the mirVANA miRNA Isolation Kit (Anibion, 
Austin, TX) and the concentration was measured spectrophoto- 
metrically (Nanodrop technologies, Montchanin, DE). The 
microRNA in the total RNA sample (80 ng/ul) was converted to 
cDNA by reverse transcription using the miRCURY LNA 
Universal RT cDNA synthesis kit (Exiqon, Vedbaek, Denmark). 
qRT-PCR using appropriate endogenous controls (SNORA66, 
UBsnRNA, RNUlAl) was performed in order to validate the 



microarray results of the top up-regulated miRNAs: miR-3648, 
miR-489, miR-638, miR-3656, miR-3687, miR-663b, miR-3126- 
5p, miR-25-5p, miR-1908, miR-204-3p, miR-711, miR-516b-5p, 
miR-498 and miR-612; as well as the top down-regulated ones: 
miR-145-3p, miR-551b-5p, miR-155-3p, miR-106b-3p, miR-140- 
5p, miR-4284, miR-520g, miR-520e, miR-944, miR-656 and 
miR-141-3p. The cDNA template was diluted 80x and amplified 
using SYBR® Green master mix (Exiqon, Vedbaek, Denmark) and 
LNA microRNA-specific primers on a ViiA 7 Real-Time PGR 
System (Applied Biosystems). Polymerase activation/DNA dena- 
turation was at 95°G for 10 mins and qPGR amplification was 
performed for 45 cycles at 95°G for 10 sec and 60°C for 1 min 
with a ramp-rate 1.6°G/sec, followed by melt curve analysis. AU 
reactions were performed in triplicates. Before proceeding with 
data normalization, a technical quality assessment was performed 
based on results of the melting curve, serial dilution curve and no- 
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Figure 3. Hierarchical Clustering (HCI). The unsupervised two-way HCI with Euclidian distance depicts differential miRNA expression in ccRCC, 
papRCC, chRCC and UT-UC. The log2 fold change in each RCC subtype and UT-UC versus the normal kidney tissue was used to construct the heat 
map. miRNA profiling accurately discriminated between RCC and UT-UC, as well as among ccRCC, papRCC and chRCC. ccRCC, clear cell renal cell 
carcinoma; papRCC, papillary renal cell carcinoma; chRCC, chromophobe renal cell carcinoma; UT-UC, upper tract urothelial carcinoma. Red and blue 
colours show significant up- or down-regulation of each miRNA in the tumour versus the normal kidney, respectively. 
doi:10.1371/journal.pone.0091646.g003 



template-controls. The stability and ranking of the endogenous 
controls was calculated with the NormFinder algorithim (www. 
mdl.dk/publicationsnormfinder.htm). For every endogenous con- 
trol gene, the pair-wise variation with all other endogenous 
controls was determined as a gene stability measurement M. 
Genes with expression stability M>1.5 were considered unstable 
across the samples and unsuitable for endogenous controls in this 
experiment. Relative expression was performed using the AACt 
method as previously described [21,27]. Data were standardized 
by log2 transformation. 

Locked nucleic acids-in situ hybridization (LNA-ISH) 

For LNA-ISH we used the probe sequences: miR-25-5p, 5'- 
DIG/CAATTGCCCAAGTCTCCGCCT/DIG-3'; miR-3687, 
5'-DIG/ACGTCGCACGAACGCCTGT/DIG-3'; LNA U6 
snRNA control probe, 5'-DIG/CACGAATTTGCGTGT- 
CATCCTT-3'; LNA scrambled microRNA control probe, 5'- 
DIG/GTGTAACACGTCTATACGCCCA/DIG-3'. Briefly, 6- 
|J,m sections of FFPE tumour and control tissue samples were 
prepared. Paraffm was melted at 60°C on the day before 
conducting the ISH experiment and slides were stored overnight 
at 4°C in an RNase-free environment. The following day, slides 
were deparaffmized in xylene and ethanol washes, treated with 
proteinase K (15 |tg/nil) for 30 minutes at 37°C, washed in PBS 
and dehydrated in ethanol. Hybridization mix (1 nM LNA U6 
snRNA, 50 nM double-DIG LNA miR-25-5p and miR-3687 
probes or 50 nM double-DIG LNA scrambled microRNA probe) 
(Exiqon, Vedbaek, Denmark) was applied for 2 hours at ~30°C 



below the RNA Tm of the corresponding LNA probe, in a 
StatSpin ThermoBrite System (Abbott Laboratories Ltd). Three 
diflFerent reactions were performed per slide. This was achieved by 
cutting with scissors the polypropylene coverslips to a size shghdy 
larger than the tissue section. Three serial sections were placed per 
glass slide, in order to perform the test, the negative and positive 
controls on the same glass slide. This not only allowed for more 
efficient work, but permitted the U6 snRNA and scrambled 
sequence controls to be done under the same exact conditions as 
miR-25-5p or miR-3687. SKdes were washed in SSC buffers prior 
to incubation with blocking solution. Subsequently, sheep anti- 
DIG-AP reagent was applied for 1 hour, slides were washed in 
PBS-T, incubated firstly with AP substrate (NBT/BCIP containing 
levamisol) for 2 hours at 30°C in a humidifying chamber protected 
from light, and then in KTBT buffer to stop the reaction. Nuclei 
were counter stained with Nuclear Fast Red. Slides were carefully 
rinsed in tap water and dehydrated in ethanol solutions. Finally, 
shdes were mounted, left overnight to allow settlement of the 
precipitate and analyzed under a light microscope (Nikon eclipse 
TE2000-U) the subsequent day. 

Computational analysis of the deregulated microRNAs 

Chromosomal location and phylogenetic analysis of DE 
miRNAs. Using miRBase vl9 (www.mirbase.org/), the chro- 
mosomal distribution of the DE miRNAs was compared with 
reported genomic alterations in the diflFerent RCC subtypes and 
UT-UC. The unaligned fasta format of the mature sequences of 
miRNA gene families was fetched using miRBase vl9 and multiple 
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Figure 4. Overlapping relationship of the deregulated miRNAs. Venn diagrams illustrate the overlapping relationship of the number of up- 
regulated mlRNAs among RCC subtypes (A), down-regulated miRNAs among RCC subtypes (B), up-regulated miRNAs among RCC subtypes and UT- 
UCs (C), down-regulated miRNAs between RCC subtypes and UT-UCs (D). Ninety-four miRNAs were co-upregulated among ccRCC, papRCC and 
chRCC; and 1 1, 44 and 24 miRNAs were specifically up-regulated in each one of the three RCC subtypes (ccRCC, chRCC and papRCC), respectively. On 
the other hand, 222 miRNAs were co-down-regulated in the three RCC subtypes, whereas 16, 18 and 5 miRNAs were specifically down-regulated in 
ccRCC, chRCC and papRCC, respectively. When the DE miRNAs in each RCC subtype were combined with those in UT-UC, we identified 89 and 206 
miRNAs that were up- and down-regulated, respectively in all tumor types. 
doi:10.1371/journal.pone.0091646.g004 



sequence alignment wa.s sub.sequently performed with the 
ClustalW algorithm using T-Coffee (www.ebi.ac.uk/Tools/msa/ 
tcofFee/). Finally, a phylogenetic tree for the miRNA sequences 
was built using bayesian reconstruction. 

miRNA target prediction. The DE microRNAs were 
imported into the miRWalk algorithm (www.ma.uni-heidelberg. 
de/ apps/ zmf/ mirwalk) and prediction of miRNA target sequences 
on the 3'UTR of their target genes wa.s performed using 7 
nucleotides as the minimum seed number. A multiple comparison 
using 4 additional algorithms was performed for filtering purposes, 
each one working based on different sets of properties regarding 
mRNA-miRNA targeting; TargetScan, miRanda, miRDB and 
Diana-microT. Filtering was performed for pairs of miRNA- 
mRNA targets that were predicted by all five algorithms and gave 
a p-value<0.05, according to the miRWalk algorithm. This p- 
value depicts the probability of a miRNA 5 '-end sequence to be 
randomly paired with a given 3'UTR mRNA sequence. The 
validated gene targets of the DE miRNAs in each RCC subtype 
and UT-UC were also investigated. 

Ingenuity pathway analysis for the significantly 
deregulated miRNAs. Ingenuity Pathway Analysis, version 7 



(IPA; Ingenuity Systems, USA; www.analysis.ingenuity.com) was 
performed to identify the molecular pathways and functional 
groupings based on published literature for the significantly DE 
miRNAs. miRNAs were uploaded into IPA and overlaid onto a 
global molecular network developed from information contained 
in the apphcation. Networks of these genes were generated by IPA 
based on their connectivity, each ranked by a score. This score is 
based on the hypergeometric distribution, calculated with the 
right-tailed Fisher's Exact Test, and corresponds to the negative 
log of this p-value. Functional analysis in IPA identified the 
published biological functions that were most significandy 
associated with the genes in the network. Genes or gene products 
are represented as nodes, where shape indicates functional groups, 
and the biological relationship between two nodes is represented as 
an edge (line). All lines are supported by at least one reference in 
literature, textbook, or from canonical information stored in the 
Ingenuity Pathways knowledge database. 

Functional assays 

Identification of anti-correlated miRNA/mRNA 
pairs. Using miRWalk, we searched for gene targets of the 
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Figure 5. ROC analysis using microarray data. ROC curves of the top 30 up-regulated miRNAs In each tissue type using the microarray 
expression data. Of them, the mIRNAs with a p<0.01 and an AUO0.8 were selected as successful distinguishing markers between cancerous and 
normal tissues. The median area under the curve (AUC) for ccRCC was 0.85 (A); for papRCC, median AUC = 0.94 (B); for chRCC, median AUC = 0.84 (C) 
and for UT-UC, median AUC = 0.94 (D). 
dol:10.1371/journal.pone.0091646.g005 



top deregulated miRNAs as described above. We also selected the 
top co-deregulated (top co-up- and top co-down-regulated) mRNA 
transcripts from 5 ccRCC GEO datasets (accession numbers: 
GSE6344, GSE4125, GSE15641, GSE781 and GSE6280). For 
functional analysis we chose the miRNA/mRNA pairs that 
exhibited anti-correlated profiles between our miRNA profiling 
and the co-deregulated mRNA transcripts from the aforemen- 
tioned 5 GEO datasets. 

Bacterial transformation of ligated products and DNA 
sequencing. Primers were designed to introduce a Spel and a 
Hindlll restriction enzyme sites to be cloned into the pMiR- 
REPORT luciferase vector (Table SI). DH5a E. coli cells 
(Takara, Japan) were gown from glycerol stocks and plasmid DNA 
was isolated and linearized using Spel and a Hindlll restriction 
enzymes. PCR-generated inserts were gel purified and ligated into 
the pMiR-REPORT luciferase vector. Ligation products were 
transformed into competent DH5a cells and the next day clones 
were picked and plasmid DNA was isolated. Insert verification 
included a restriction reaction with Spel and Hindlll. Also, 1 00 ng 
of plasmid DNA containing the insert was sequenced to confirm 
the absence of unwanted mutations, using the forward sequencing 
primer: 5'-AGGCGATTAAGTTGGGTA-3'. 

DNA sequencing of predicted target regions was performed 
using BigDye V3.1 chemistiy on an ABI Prism Genetic Analyzer 
(Applied Biosystems, OA). Sequencing primers were designed to 
flank the target region but also included an additional 300 bp on 
average on each side. Sequence electropherograms were obtained 
from the ABI Sequencing Analysis V5.2 software (Applied 
Biosystems, CA) and sequences were imported into MEGA4 
software to be aligned against a reference sequence with ClustalW 
algorithm. 



Real-Time PCR for endogenous mlRNA detection 

To examine the expression levels of the endogenous miR-23c, 
mR-204-3p, miR-520a-5p, miR-489, miR-874 and miR-513a-3p 
in undifferentiated AB8/13 cells, total RNA enriched in small 
RNAs was isolated from the cells using TRIzol® reagent (Ambion). 
miRNA specific reverse transcription was performed with the 
miRCURY LNA Universal RT cDNA synthesis kit (Exiqon, 
Vedbaek, Denmark). The cDNA template was diluted 80x and 
amplified using SYBR® Green master mix (Exiqon, Vedbaek, 
Denmark) and LNA microRNA-specific primers on a ViiA 7 Real- 
Time PCR System (Applied Biosystems). AH reactions were 
performed in triphcates. SNORA66, U6 snRNA and RNUlAl 
were used as reference. 

Transfection of AB8/13 podocytes and luciferase assays 

To evaluate the binding efficiency of the top deregulated 
miRNAs onto their predicted targets, the pMiR-REPORT miRNA 
Expression Reporter Vector System (AMBION, Texas) was used. 
The AB8/13 undifferentiated podocyte cells were incubated at 
33°C at 5% CO2 and cultured in RPMI medium, supplemented 
with 10% Fetal Bovine Serum (FBS) (Invitrogen, CA), 1% of 
100 units/ml Penicillin/ Streptomycin (Invitrogen, CA) and 1% 
Insulin-Transferrin-Selenium (Invitrogen, CA). For the luciferase 
reporter system experiments, AB8/13 cells were transfected with 
equal amounts of the pMIR- REPORT luciferase and P-gal vectors 
(100 ng each), as well as with 25 nM of miRNA inhibitors, mimics, 
or the AUStars Negative Control scrambled sequence LNA (Qjagen, 
West Sussex, UK), using lipofectamine 2000 (Invitrogen, CA). The 
P-gal vector was used for normalization. Every experiment was 
performed in quadruplicates in 12-weIl cell culture plates with the 
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Table 2. ROC statistics of the mlRNAs that exhibited the best discriminatory ability in ccRCC (A), papRCC (B), chRCC (C) and UT-UC 
(D) vs. the normal kidney tissue. 



A. miRNAs in ccRCC 


AUC 


S.E. 


95% 


C.i. 


p value 


Sensitivity 


1 -Specificity 


miR-3648 


0.918 


0.057 


0.746- 


0.971 


0.0001 


100.0 


90.0 


miR-143-3p 


0.878 


0.060 


0.735- 


0.969 


0.0001 


100.0 


95.0 


miR-25-5p 


0.878 


0.066 


0.723- 


0.982 


0.0001 


100.0 


95.0 


miR-628-3p 


0.868 


0.064 


0.724- 


0.975 


0.0001 


100.0 


90.0 


miR-921 


0.859 


0.060 


0.727- 


0.964 


0.0002 


100.0 


95.0 


miR-210 


0.853 


0.066 


0.699- 


0.958 


0.0004 


100.0 


95.0 


miR-711 


0.853 


0.071 


0.678- 


0.958 


0.0006 


100.0 


95.0 


miRPIus-C1087 


0.850 


0.067 


0.684- 


0.947 


0.0007 


100.0 


94.7 


miR-516b-5p 


0.846 


0.069 


0.674- 


0.944 


0.0009 


100.0 


95.0 


miR-3126-5p 


0.828 


0.072 


0.660- 


0.944 


0.0012 


100.0 


90.0 


miR-3687 


0.818 


0.072 


0.661- 


0.941 


0.0013 


100.0 


95.0 


miR-204-3p 


0.816 


0.075 


0.650- 


0.946 


0.0014 


100.0 


90.0 


miR-3676-3p 


0.809 


0.074 


0.650- 


0.941 


0.0015 


100.0 


90.0 


miR-489 


0.802 


0.075 


0.644- 


0.939 


0.0018 


100.0 


95.0 


miR-498 


0.801 


0.072 


0.645- 


0.928 


0.0022 


100.0 


90.0 


B. mIRNAs in papRCC 


AUC 


S.E. 


95% 


CI. 


p value 


Sensitivity 


1 -Specificity 


miR-4290 


0.985 


0.035 


0.901- 


1.042 


0.002 


100 


92.8 


miR-3687 


0.971 


0.057 


0.829- 


1.056 


0.004 


100 


92.8 


miRPIus-C1087 


0.957 


0.052 


0.840- 


1.046 


0.004 


100 


92.8 


miR-3648 


0.942 


0.063 


0.803- 


1.054 


0.005 


100 


92.8 


miR-3676-3p 


0.928 


0.043 


0.871- 


1.043 


0.003 


100 


92.8 


miR-3195 


0.928 


0.022 


0.941- 


1.030 


0.001 


100 


92.8 


miR-630 


0.871 


0.035 


0.901- 


1.042 


0.002 


100 


92.8 


miR-3126-5p 


0.857 


0.059 


0.812- 


1.045 


0.005 


100 


92.8 


C. miRNAs in chRCC 


AUC 


S.E. 


95% 


CI. 


p value 


Sensitivity 


1 -Specificity 


miR-1276 


0.944 


0.122 


0.516- 


0.994 


0.1731 


100 


93.33 


miRPIus-C1087 


0.911 


0.133 


0.538- 


1.061 


0.109 


100 


93.33 


miR-3665 


0.911 


0.138 


0.573- 


1.116 


0.066 


100 


93.33 


miR-1909-3p 


0.889 


0.133 


0.538- 


1.061 


0.109 


100 


93.33 


miR-204-3p 


0.867 


0.138 


0.485- 


1.026 


0.173 


100 


93.33 


miR-193b-5p 


0.844 


0.179 


0.381- 


1.085 


0.213 


100 


93.33 


miR-3648 


0.822 


0.169 


0.356- 


1.021 


0.314 


100 


93.33 


miR-4290 


0.822 


0.163 


0.414- 


1.052 


0.213 


100 


93.33 


miRPIus-C1076 


0.822 


0.074 


0.764- 


1.057 


0.028 


100 


93.33 


miR-663b 


0.800 


0.198 


0.255- 


1.034 


0.441 


100 


93.33 


miR-489 


0.800 


0.194 


0.329- 


1.092 


0.26 


100 


93.33 


D. miRNAs in UT-UC 


AUC 


S.E. 


95% 


CI. 


p value 


Sensitivity 


1 -Specificity 


miR-3648 


1.00 


0.063 


0.818- 


1.066 


0.009 


100 


92.308 


miR-489 


1.00 


0.078 


0.768- 


1.078 


0.001 


100 


92.308 


miR-3665 


1.00 


0 


1 .000- 


1.000 


0.003 


100 


92.308 


miR-874 


1.00 


0.063 


0.818- 


1.066 


0.009 


100 


92.308 


miR-4285 


1.00 


0 


1 .000- 


1.000 


0.003 


100 


92.308 


miR-210 


1.00 


0 


1 .000- 


1.000 


0.003 


100 


92.308 


miR-612 


0.98 


0.068 


0.788- 


1.058 


0.012 


100 


92.308 


miR-663b 


0.96 


0.046 


0.869- 


1.054 


0.006 


100 


92.308 


miR-3687 


0.96 


0.045 


0.872- 


1.051 


0.006 


100 


92.308 


miR-1276 


0.96 


0 


1 .000- 


1.000 


0.003 


100 


92.308 


miR-1908 


0.94 


0.056 


0.831- 


1.053 


0.009 


100 


92.308 
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miR-3195 


0.94 


0.087 


0.71 3-1 .055 


0.023 


1 00 


92.308 


miRPIus-C1087 


0.94 


0 


1.000-1.000 


0.003 


100 


92.308 


miR-498 


0.94 


0.063 


0.81 8-1 .056 


0.009 


1 00 


92.308 


miR-14fiQ 


0 94 


0 091 


n fiRfi-1 044 


0 031 


1 00 


92 308 


1 1 1 1 rv 


0 94 


0 094 


n 71 8-1 ORQ 


0 017 


1 00 


92 308 


[ 1 1 1 rv 3 1 _7\J 


0 92 


0 124 


n 607-1 090 


0 041 


1 00 


92 308 


miR-204-3p 


0.92 


0.045 


0.872-1 .051 


0.006 


1 00 


92.308 


miR-1 85-3p 


0.91 


0.06 


0.824-1 .060 


0.009 


1 00 


92.308 


miR-423-5p 


0.90 


0.1 22 


0.567-1 .048 


0.07 


1 00 


92.308 


miR-921 


0.88 


0.087 


0.712-1.057 


0.023 


1 00 


92.308 


miR-638 


0.86 


0 


1 .000-1 .000 


0.003 


1 00 


92.308 


miR-25-5p 


0.86 


0.029 


0.922-1.039 


0.004 


100 


92.308 


miR-3676-3p 


0.84 


0.071 


0.774-1.053 


0.014 


100 


92.308 


miR-71 1 


0.84 


0.046 


0.869-1.054 


0.006 


100 


92.308 


miR-3656 


0.82 


0 


1.000-1.000 


0.003 


100 


92.308 


miR-516b-5p 


0.80 


0.063 


0.818-1.066 


0.009 


100 


92.308 





Due to the small chRCC sample number (n = 3), no statistical significance could be gained. AUC, area under the curve; S.E., standard error; 95% C.I., 95% confidence 
interval. 

doi:l 0.1 371 /journal.pone.0091 646.t002 



appropriate controls. Cells were harvested 16 hours after transfec- 
tion. The Dual-Light Assay Kit (Applied Biosystems, CA) was used 
for the quantification of both luciferase and P-gal in an automated 
luminometer (Sirius, Berthold Detection Systems, Pforzheim, 
Germany). 

Western blotting 

AB8/13 cells were lysed in equal volumes of pre-heated 2xSDS 
loading buffer (Sodium Dodecyl Sulphate-125 mM Tris-HCl 
pH 6.8, 20% Glycerol, 2% SDS, 2% P-mercaptoethanol and 
bromophenol blue) and homogenized using a 2 ml syringe. Whole 
cell lysates were subsequently electrophoresed in a 12% SDS- 
Polyacrylamide gel. Gel transfer was held in a wet transfer system 
on Hybond Polyvinylidene Fluoride (PVDF-Millipore, MA) 
membranes. Membranes were blocked with 5% non-fat dry milk 
in PBS/0.01% Tween20 for 1 hour at room temperature. Primary 
antibody was diluted in milk and added to the membrane for one 
hour. TFCP2L1 protein was detected with the goat primary 
polyclonal antibody C-20 (SantaCruz Biotechnology, CA) at 
around 48 kDa. P-Tubulin was used as loading control by using 
the T-4026 primary antibody (Sigma, Taufkirchen, Germany). As 
secondary antibody we used the rabbit anti-mouse antibody 
(SantaCruz Biotechnology, CA), conjugated with Horseradish 
Peroxidase (HRP). Proteins were detected using the Enhanced 
ChemiLuminescence (ECL) Plus Blotting Detection .system 
(Amersham Biosciences, Buckinghamshire, UK) and were visual- 
ized on a ChemiDoc MP imaging system using the Image Lab 
software (Bio-Rad Laboratories). 

Statistical analyses 

Microarray data analysis was performed using Matlab®. The 
two tailed student t-test was used to test the mean differences 
between two groups. Continuous variables are expressed as 
median ± standard deviation unless indicated differently. Receiver 
operating characteristic (ROC) curves were established to evaluate 
the diagnostic value of deregulated miRNAs for differentiating 
between RCC subtypes and controls. Pearson's correlation 
coefficient (CC) was used to identify positive or negative linear 
correlation between variables. Linear regression analysis using 
miRNA expression levels as dependent variables and patient age. 



sex, tumour stage, grade, size and nodal status as independent 
variables, was also performed to identify miRNAs associated with 
the aforementioned characteristics. Linear regression analysis 
along with ANOVA were also performed to develop models that 
discriminate RCC/UT-UC from the normal kidney samples. The 
median expression values of all the DE miRNAs were used to 
perform linear regression analysis. Similar analysis was performed 
based on tumor grade. All tests were performed 2-sided and a 
p^0.05 was considered as statistically significant. Luciferase 
expression levels were analyzed using one-way ANOVA followed 
by Tukey post-testing, after being normalized against P-gal 
expression levels. One-way ANOVA was also used to test 
densitometry results from western blot analysis, followed by Tukey 
post-testing. All statistical analyses were performed using Graph- 
Pad Prism V.5 and Matlab®. 

Results 

microRNA profiling 

microRNA profiling was performed retrospectively in a total of 
27 RCCs (3 chRCC, 5 papRCC and 19 ccRCC), 4 UT-UCs and 20 
normal kidney tissues. We identified 434 deregulated (DE) miRNAs 
(p<0.05 and FDR<0.05) in all RCC subtypes and UT-UC versus 
the normal kidney. The majority of them was down-regulated in 
RCCs but over-expressed in UT-UCs. OveraU, 126 miRNAs (29%) 
showed increased expression and 303 miRNAs (69.8%) had 
decreased expression in RCCs and UT-UCs vs. the normal kidneys 
(Figure 2). Specifically, 374, 420, 421 and 409 miRNAs were 90% 
consistently down-regulated in ccRCC, papRCC, chRCC and UT- 
UC, respectively. Unsupervised two-way hierarchical clustering 
(HCl) with Euclidian distance accurately discriminated between 
RCC and UT-UC. Apart from one chRCC sample that was 
clustered with ccRCCs, HCl also managed to successfully classify 
the 3 RCC subtypes among them. Furthermore, it showed that 
ccRCC is more closely related to papRCC and that both are distinct 
from chRCC or UT-UC (Figure 3). HCl pair-tissue-wise also 
discriminated perfectiy ccRCC from papRCC, chRCC from 
ccRCC, chRCC from papRCC, as weU as ccRCC from UT-UC, 
papRCC from UT-UC and chRCC from UT-UC (Figure SI). 
Table S2 lists the top deregulated miRNAs for each RCC subtype 
and UT-UC. The overlapping relationship of the deregulated 
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Figure 6. qRT-PCR validation of tfie top 25 deregulated miRNAs. The Volcano-plots depict the significantly deregulated miRNAs in ccRCC, 
papRCC, chRCC and UT-UC. Eight miRNAs were significantly up-regulated in ccRCC, 5 in papRCC, 3 in chRCC and 4 in UT-UC. On the other hand, miR- 
656 was significantly down-regulated in papRCC; miR-155-3p, miR-106b-3p, miR-140-5p and miR-656 were significantly down-regulated in chRCC; 
and miR-520g was significantly down-regulated in UT-UC. The threshold of statistically significant difference was set at p<0.05 and log2 fold 
change>2. 

doi:10.1371/journal.pone.0091646.g006 



miRNAs among the three RCC subtypes, as well as among all 
tumour types is depicted in Figure 4. Tissue-specific miRNAs, 
including those that can distinguish RCCs from UT-UC, were 
identified. We found 16 miRNAs that were specifically down- 
regulated (miR-3144-5p, miR-193b-3p, miR-587, miR-3117-3p, 
miR-769-5p, niiR-617, miR-3615, miR-148b-5p, miR-24-2-5p, 
miR-142-5p, miR-3164, miR-514b-5p, miR-375, miR-520c-3p, 
miR-191-5p and miR-520b/miR-520c-3p) and 5 miRNAs that 
were specifically up-regulated in UT-UC (miR-885-5p, miR-1912, 
miR-378a-5p, miR-24-3p and miR-26a-2-3p) Table S3. 

ROC analysis of the top 30 up-regulated miRNAs in each tissue 
type, using the microarray expression data, evaluated the extent to 
which they could separate each tumour entity from the normal 



kidney tissue. Of them, the miRNAs with a p<0.01 and an 
AUO0.84 were selected as successful distinguishing markers 
between each tumor subtype and the normal kidney tissue 
(Figure 5). The ROC statistics of the miRNAs that exhibited 
the best discriminatory ability are depicted in Table 2. 

Validation of deregulated miRNA expression by qRT-PCR 

We verified by qRT-PCR the expression levels of the top 25 
deregulated miRNAs as acquired by miRNA profiling (14 most 
up-regulated and 1 1 most down-regulated miRNAs). These 
included: miR-3648, miR-489, niiR-638, miR-3656, miR-3687, 
miR-663b, miR-3126-5p, miR-25-5p, miR-1908, miR-204-3p, 
miR-71 1, niiR-516b-5p, niiR-498 and miR-612 (top up-regulated 
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Figure 7. ROC analysis using qRT-PCR data. ROC curves of the significantly deregulated mIRNAs, using the expression data from qRT-PCR 
analysis. Of them, the mIRNAs with a p<0.01 and an AUO0.8 were selected as successful distinguishing markers between cancerous and normal 
kidney tissues. The median area under the curve (AUC) for ccRCC was 0.802 (A); for papRCC, median AUC = 0.756 (B); for chRCC, median AUC = 0.926 
(C) and for UT-UC, median AUC = 0.955 (D). 
doi:1 0.1 371/journal.pone.0091 646.g007 



miRNAs); as well as miR-145-3p, miR-551b-5p, miR-155-3p, 
miR-106b-3p, miR-140-5p, miR-4284, miR-520g, miR-520e, 
miR-944, miR-656 and miR-141-3p (top down-regulated miR- 
NAs). The validated significantly deregulated miRNAs in ccRCC, 
papRCC, chRCC and UT-UC are depicted in Figure 6. 

The frequencies of differential expression were mostly consistent 
across the same subtypes. In keeping with the microarray results, 
some miRNAs showed a similar expression pattern between the 
subtypes. For example, miR-3687 exhibited a similar up-regula- 
tion pattern in all RCC subtypes and in UT-UC; miR-3656, miR- 
638 and miR-3648 were similarly up-regulated in ccRCC, 
papRCC and UT-UC and miR-25-5p was similarly up-regulated 
in ccRCC and papRCC. Also, miR-656 was equally down- 
regulated in papRCC and chRCC. 

In order to further validate the capability of the DE miRNAs to 
discriminate between tumour and normal kidney tissue, ROC 
curves were further constructed for the deregulated miRNAs in 
each RCC subtype and in UT-UC, using the normalized qRT- 
PCR data. miRNAs with a p<0.01 and an AUO0.8 were then 
selected as successful distinguishing markers between cancerous 
and normal kidney tissue (Figure 7). miR-3648, miR-638, miR- 
3656 and miR-3687 were excellent discriminators between 
ccRCC and normal tissue (median AUO0.8, p<0.0001); miR- 
3687 and miR-25-5p were very good discriminators between 
papRCC and normal tissue (median AUO0.8, p<0.05); mlR- 
3687, miR-4284 and miR-141-3p were very good discriminators 



between chRCC and normal tissue (median AUO0.8, p<0.05); 
and miR-3648, miR-3656 and miR-3687 discriminated perfectiy 
UT-UC from the normal tissue (median AUO0.94, p<0.001). 
ROC curve data are shown in Table 3. The qRT-PCR results 
revealed similar expression patterns with the microarray results 
(Figure 8). 

A blinded independent validation set composed of 40 freshly 
frozen ccRCC, papRCC, chRCC, UT-UC and normal kidney 
samples was also used to vaUdate the discriminatory ability of the 
qRT-PCR verified deregulated miRNAs. The results were highly 
repeatable between the two cohorts, recapitulating the specificity 
of these miRNAs in each tumour group (in ccRCC, Pearson's 
CC = 0.987, p<0.001; in papRCC, Pearson's CC = 0.917, 
p = 0.001; in chRCC, Pearson's CC = 0.952, p<0.001; in UT- 
UC, Pearson's CC = 0.926, p = 0.001) (Figure 9). 

miRNA expression and clinical characteristics of the 
patients 

Only a few correlations between the expression levels of the DE 
miRNAs (verified by qRT-PCR) with the tumour stage, grade and 
size were identified. The expression levels of the significantiy 
deregulated miRNAs that were verified by qRT-PCR, were 
correlated with the tumour stage, grade and size. Grade 4 ccRCCs 
exhibited higher miR-3656 and miR-3687 levels vs. those of grade 
3 (p = 0.028 and p = 0.0497, respectively; unpaired t-test). 
Furthermore, grade 4 ccRCCs exhibited higher mrR-21-5p levels 
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Table 3. ROC statistics of the miRNAs that exhibited the best discriminatory ability in ccRCC (A), papRCC (B), chRCC (C), and UT-UC 
(D) vs. the normal kidney tissue, respectively. 





A. miRNAs in ccRCC 


AUC 


S.E. 


95% CI. 


p value 


Sensitivity 


1 -Specificity 



miR-3648 0.839 

miR-489 0.773 

miR-638 0.807 

miR-3656 0.817 

miR-3687 0.847 

miR-25-5p 0.796 

miR-21-5p 0.737 

miR-663b 0.668 

B. mIRNAs in papRCC AUC 
miR-3648 0.745 
miR-638 0.650 
miR-3656 0.725 
miR-3687 0.800 
miR-25-5p 0.865 
miR-656 0.770 

C. miRNAs in chRCC AUC 
miR-3687 0.933 
miR-4284 0.891 
miR-141-3p 0.925 

D. miRNAs in UT-UC AUC 
miR-3648 0.940 
miR-638 0.905 
miR-3656 0.970 
miR-3687 0.995 



0.049 

0.086 

0.056 

0.055 

0.048 

0.054 

0.065 

0.068 

S.E. 

0.108 

0.111 

0.110 

0.094 

0.102 

0.174 

S.E. 

0.053 

0.071 

0.084 

S.E. 

0.041 

0.052 

0.027 

0.007 



0.742-0.935 

0.602-0.943 

0.697-0.917 

0.707-0.926 

0.752-0.941 

0.689-0.902 

0.608-0.866 

0.533-0.802 

95% CI. 

0.532-0.957 

0.431-0.868 

0.508-0.942 

0.614-0.985 

0.663-1.066 

0.427-1.113 

95% CI. 

0.828-1.038 

0.752-1.031 

0.759-1.092 

95% CI. 

0.857-1.022 

0.801-1.008 

0.916-1.024 

0.979-1.011 



<0.0001 

0.003 

<0.0001 

<0.0001 

<C.C001 

<0.0001 

0.001 

0.023 

p value 

0.076 

0.278 

0.104 

0.030 

0.008 

0.058 

p value 

0.013 

0.025 

0.033 

p value 

0.001 

0.003 

0.0006 

0.0003 



100.0 
100.0 
100.0 
100.0 
100.0 
100.0 
100.0 
100.0 

Sensitivity 

100 
100 
100 
100 
100 
100 

Sensitivity 

100 
100 
100 

Sensitivity 

100 
100 
100 
100 



97.5 

96.66 

97.5 

97.5 

97.5 

97.5 

96.0 

97.5 

1 -Specificity 

97.5 
97.5 
97.5 
97.5 
97.5 
97.5 

1 -Specificity 

97.5 
97.5 

88.89 

1 -Specificity 

97.5 
97.5 
97.5 
97.5 



AUC, area under the curve; S.E., standard error; 95% C.I., 95% confidence interval. 
doi:1 0.1 371 /journal.pone.0091 646.t003 



VS. grade 2 ccRCCs (p = 0.0449, unpaired t-test). We also found a 
significant association between miR-21-5p expression and ccRCC 
size (Pearson's CC = 0.442, p = 0.031). The levels of miR-25-5p 
were significantly associated with papRCC tumour size (Pearson's 
CC = 0.935, p = 0.019). Grade II UT-UCs expressed higher levels 
of miR-3648 and miR-638 vs. grade III ones (p = 0.0153 and 
p = 0.0057, respectively). Finally, niiR-638 was correlated with 
tumour size (Pearson's CC = 0.878, p = 0.050). 

Furthermore, miRNA expression was validated using one-way 
ANOVA along with linear regression. Of note, miR-204-3p, miR- 
25-5p, miR-3648, miR-3687, miR-489, miR-516b-5p, miR-638 
and miR-663b exhibited a linear expression pattern among the 3 
RCC subtypes, UT-UC and the normal kidney samples {R~ 
range = 0.93-0.98). Also, linear regression could classify between 
RCC and UT-UC, with the latter expressing higher levels of the 
pre-mentioned miRNAs (Figure S2). 

Similarly, high linearity was achieved when we grouped 
together ccRCC, papRCC and chRCC samples, and we 
compared them to the UT-UC and the normal kidney cases 
(miR-204-3p, R'>0.99; miR-25-5p, if^>0.83; miR-3648, 
i^'>0.98; miR-3687, i?^>0.95; miR-489, i^'>0.99; and miR- 
638, if^>0.98) (Figure S3). 

Linear regression of RCC grades with respect to miRNA 
expression, revealed that miR-148b-5p, miR-1910 and miR-200b- 
5p, miR-210 and miR-3656, presented linear behavior in a 
descending order of the RCC grade; whereas miR-3687 and miR- 



654-5p manifested linear behavior in ascending order of the RCC 
grade (Figure S4). 

Furthermore, miR-148a-5p and miR-940 manifested excellent 
linear behavior in an descending order of age {R^ — 0.96 and 0.99, 
respectively), pointing out that both might be age-specific. Also, 
regression analysis showed that miR-145-3p and miR-1973 
manifested increasing expression levels with ascending tumor size 
(<5cm, 5-lOcm and 10-15 cm) (i?- = 0.95 and 0.92, respectively). 

Other varieties, such as a patient's sex, age or tumor stage, were 
not associated with any other miRNA expression levels. 

Pair-wise miRNA correlations 

Several DE microRNAs were correlated among each other. In 
ccRCC, the statistically significant Spearman's correlation coeffi- 
cients ranged from rs = -0.585 (miR-106b-3p to miR-612; 
p = 0.002) to r, = 0.843 (miR-638 to miR-3687; p<0.0001). In 
papRCC, the statistically significant Spearman's correlation 
coefficients ranged from rs = -0.90 (miR-145-3p to mLR-711; 
and miR-21-5p to miR-141-3p; both p = 0.083) to r,= 1.0 (miR- 
638 to mLR-3648; miR-663b to miR-3656; miR-516b-5p to miR- 
489; miR-498 to miR-489; miR-155-3p to miR-489, among 
others; all p = 0.017). Furthermore, in UT-UC the statistically 
significant Spearman's correlation coefficients ranged from r^ — 
-0.9 (miR-551b-5p to miR-3126-5p; miR-520g to miR-3126-5p; 
miR-141-3p to niiR-1908, among others; aU p = 0.083) to r,= 1.0 
(miR-489 to miR-3648; miR-3656 to miR-3648; miR-3656 to 
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Figure 8. Correlation between microarrays and qRT-PCR. Median log2 fold change expression levels of the 14 most up-regulated and 1 1 most 
down-regulated mlRNAs between ccRCC, papRCC, chRCC and UT-UC and the normal kidney tissue, as determined by both qRT-PCR and microarray 
analysis. As shown in the figure, qRT-PCR and microarray results were highly compatible. The most identical results between the two techniques were 
those for ccRCC, which was expected due to the high sample number (Pearson's CC = 0.778, p<0.001). The qRT-PCR results for papRCC, chRCC and 
UT-UC also revealed similar deregulation patterns with those of the microarray experimentation, however the correlation coefficients were lower, 
apparently due to small sample number (in papRCC, CC = 0.596, p = 0.002; in chRCC CC = 0.570, p = 0.003; in UT-UC, CC = 0.517, p = 0.009). 
doi:10.1371/journal.pone.0091646.g008 



miR-489; miR-155-3p to miR-638; miR-122-5p to miR-638, 
among others; all p = 0.017). Finally, in the normal kidney the 
statistically significant Spearman's correlation coefficients ranged 
from r, = -0.325 (miR-551b-5p to miR-155-3p; p = 0.041) to 
r, = 0.885 (miR-3648 to miR-3656; p<0.0001). 

Validation of the pattern of expression of miR-25-5p and 
miR-3687 by LNA-ISH 

In situ hybridization using digoxigenin labeled LNA probes 
permitted the spatial localization of two of the most deregulated 
miRNAs: a) mLR-25-5p which was significantly up-regulated in 
ccRCC and papRCC; and b) miR-3687 which was significantly 
up-regulated in all 3 RCC subtypes as well as in UT-UC. The 
detection of both miRNAs was confmed to the cytoplasm. The 
nuclear expression of U6 snRNA was confirmed in all patient 
samples, whereas the scrambled oligonucleotide was negative in all 
samples. Each RCC or UT-UC section was compared against its 
corresponding normal kidney section. The elevated expression of 
miR-25-5p was confirmed in all ccRCC sections. Furthermore, 



ccRCCs of high stage and grade stained stronger miR-25-5p 
compared to lower stage and grade tumours. Also, papRCCs of 
type II stained stronger for miR-25-5p compared to type I 
papRCCs (Figure 10). Likewise, the elevated expression of miR- 
3687 was confirmed in the majority of both RCC and UT-UC 
sections. 

Correlation of miRNA deregulation and chromosomal 
aberrations in RCC and UT-UC 

The chromosomal locations of the DE miRNAs were compared 
to the most commonly reported chromosomal aberrations in 
ccRCC, chRCC, papRCC [28] or UT-UC [29]. In general, 
chromosomes 14q, 19q, Xq exhibited the highest number of 
deregulated miRNAs. In each RCC subtype, miRNAs had 
deregulated patterns that agreed with some of the previously 
reported chromosomal gains and losses; i.e., up-regulated miRNAs 
were located in regions of chromosomal gains and vice versa. 
Table S4 lists a summary of the reported chromosomal 
abnormalities associated with the three RCC subtypes and UT- 
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Figure 9. Validation in a blinded independent cohort. A blinded independent validation cohort composed of 40 freshly frozen ccRCC, papRCC, 
chRCC, UT-UC and normal kidney samples was also used to validate the discriminatory ability of the qRT-PCR verified deregulated miRNAs. The results 
were highly repeatable between the two cohorts, recapitulating the specificity of these miRNAs in each tumour group. 
doi:10.1371/journal.pone.0091646.g009 



UC, as well as some of the well-studied genes located within those 
aberrated regions that are involved in the pathogenesis of the 
tumour subtypes. Many of the significantly deregulated miRNAs 
are located within cytogenetic regions that are significantly 
aberrated in the different subtypes. Chromosomal mapping for 
the deregulated miRNAs in each RCC subtype and in UT-UC is 
depicted in Figure 11. 

Phylogenetic analysis of the DE miRNAs in RCC 

Phylogenetic analysis showed that 340 of the 434 (78.3%) DE 
miRNAs were highly consei"ved among species. Specifically, 
70.68% of the up-regulated miRNAs and 81.22% of the down- 
regulated miRNAs were phylogeneticaUy conserved among 
species, respectively. Figure S5 shows representative multiple 
sequence alignments and trees of miR-489 and miR-25-5p. 

IPA Core Analysis was performed for the significantly DE 
miRNAs, in each RCC subtype and in UT-UC. The loga ratio 
was set at 0.3 (range -1.9155 to 1.1864) and the DE miRNAs were 
filtered down to 61 for ccRCC, 48 for papRCC, 44 for chRCC 



and 77 for UT-UC. The top biological functions and the major 
miRNA networks are shown in Figure 12. 

Target prediction and functional analysis of anti- 
correlated mRNA/miRNA pairs 

From the GEO datasets with accession numbers GSE6344, 
GSE4125, GSE 15641, GSE781 and GSE6280 we isolated the top 
co-up- and down-regulated genes in ccRCC. We also isolated the 
top co-upregulated and top co-downregulated microRNAs 
according to our microarray analysis in ccRCC. Using the 
miRWalk algorithm (http://www.ma.uni-heidelberg.de/apps/ 
zmf/mirwalk) and four other prediction algorithms (TargetScan, 
miRanda, miRDB and Diana-microT) we identified the top anti- 
correlated mRNA/miRNA pairs, with a seed length (SL) of 10 
nucleotides or more and a p<0.001. Using this kind of analysis we 
concluded to the following anti-correlated mRNA/ miRNA pairs: 
ABCAl/miR-23c (SL=10; p = 0.0032), ATP2B4/miR-520a-5p 
(SL=10; p = 0.004) and NR3Cl/miR-513a-3p (SL=11; 
p = 0.0009) (up-regulated mRNA/ down-regulated miRNA); and 
TFCP2L1 /miR-489 (SL=10; p = 0.0074), TRIMlO/miR-874 
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Figure 10. Locked Nucleic Acids-ln Situ Hybridization (LNA-ISH). LNA-ISH for miR-25-5p in ccRCC (A), papRCC (B), chRCC (C) and UT-UC (D). 
miR-25-5p was confined to the cytoplasm both in normal and tumour sections. The nuclear expression of U6 snRNA (positive control) was confirmed 
in all the patient samples, whereas the scrambled oligonucleotide was negative in all samples. miR-25-5p high expression was confirmed in all ccRCC 
sections by LNA-ISH. ccRCCs of high stage and grade stained stronger miR-25-5p compared to lower stage and grade ccRCCs. Each ccRCC section was 
compared against its corresponding normal kidney section (A). papRCCs of type II stained stronger for miR-25-5p compared to type I papRCCs (B). 
Validating the qRT-PCR results, miR-25-5p did not stain stronger in chRCC sections vs. the normal tissue ones. This was also confirmed for chRCCs 
with focal sarcomatoid differentiation, suggesting that miR-25-5p does not play any role in the metastatic behavior of the tumour (C). Verifying the 
qRT-PCR results, miR-25-5p was not significantly stronger in UT-UC vs. the normal tissue sections (D). 
doi:10.1371/journal.pone.0091646.g010 



{SL=10; p = 0.002), SLCA7A8/miR-210 (SL=10; p = 0.004) 
(down-regulated mRNA/up-regulated miRNA). 

The AB8/13 undifferentiated podocyte cells were incubated at 
33°C at 5% C02 and cultured in RPMI medium, supplemented 
with 10% Fetal Bovine Serum (FBS), 1 % of 100 units/ml Penicillin/ 
Streptomycin and 1 % Insuliti-Transferrin-Selenium (ITS). Endog- 
enous expression of miR-23c, miR-520a-5p, miR-513a-3p, miR- 
489, miR-874 and miR-210 in undifferentiated AB8/13 podocyte 
cells was verified by qRT-PCR. A segment of the 3'UTR of the top 
co-deregulated genes was cloned into the 3'UTR of the luciferase 



gene in pMIR-REPORT plasmid (pMIR-REPORT-Gene of 
interest). Luciferase reporter plasmids and P-gal reference plasmid 
were co-transfected in AB8/13 podocyte cell line with the 
corresponding miRNA mimic, miRNA inhibitor or negative control 
mimics for 1 6 hours followed by luciferase and P-gal measurement. 
Co-transfection of pMIR-REPORT-TFCP2L 1 with either miR- 
489 mimic or inhibitor resulted in significant reduction or increase 
in luciferase expression versus the negative control, respectively, 
indicating that miR-489 binds directly on the 3'UTR of TFCP2L1 
(p<0.001; ANOVA). In agreement, transfection of miR-489 mimics 
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Figure 11. Correlation of miRNA deregulation and chromosomal aberrations. Chromosomal mapping for the deregulated mlRNAs in 
ccRCC (A), papRCC (B), chRCC (C) and UT-UC (D). In ccRCC, the up-regulated miRNAs were mainly mapped on chromosomes Xq, 17q, 5q, 14q, 19q, 
11q, 12q, 3p and lOp (in descending order); whereas the down-regulated ones were mainly mapped on chromosomes 14q, 19q, Xq, 8q, 1p, 17q, 7q, 
3q, lOq, 12q and 2q (A). In papRCC, the up-regulated miRNAs were mainly mapped on chromosomes 19q, 14q, Xq, 5q, 3p, 16p, llq, 9q and 15q; 
whereas the down-regulated miRNAs were mainly mapped on chromosomes 14q, 19q, Xq, 1 p, 7q, 8q, 12q, lOq, 17q, 3q, 19p, Xp, 2q, 5, 6q, 4q, 1 1q, 
15q and 20q (B). In chRCC, the up-regulated miRNAs were mainly mapped on chromosomes 14q, Xq, 19q, 17q, 5q, 7q, 1 1q, 15q, 2q and 12q; whereas 
the down-regulated ones were mainly mapped on chromosomes 1 4q, 1 9q, Xp, 3p, 8q, 1 7q, 3q, 1 0q, 1 2q, 4q, 7q, 3p, 5q, 1 9p and 20q (C). In UT-UC, the 
up-regulated miRNAs were mostly mapped on chromosomes 1 7q, 1 9q, llq, 1 6p, 5q, 1 4q, 1 5q, Xq, 1 p and 3p; whereas the down-regulated miRNAs 
were mainly mapped on chromosomes 14q, 19q, Xq, 8q, 1p, 7q, 12q, 3q, lOq, 5q, 17p, 19p, Xp, 2q, 4q, 3p and 8p (D). 
doi:10.1371/journal.pone.0091646.g011 
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Figure 12. Ingenuity pathway analysis (IPA). A. In ccRCC, tine most 
important biological functions of the top deregulated mlRNAs were: 1) 
Cancer (p = 8.71 E-1 1-4.91 E-02); 2) Renal and Urological Disease 
(p = 8.71 E-11 -2.1 2E-09); 3) Inflammatory Disease (p = 2.12E-09-2.14E- 
03); 4) Inflammatory Response (p = 2.12E-09-2.4E-04) and 5) Reproduc- 
tive System Disease (p = 3.37E-05-2.58E-02). The associated functions of 
the major miRNA network (score = 30) were: Inflammatory disease, 
inflammatory response, renal inflammation. Argonaute RISC catalytic 
component 2 (EIF2C2), tumour protein p53 (TP53), v-myc myelocyto- 
matosis viral oncogene homolog (MYC) and EPH receptor B6 (EPHB6) 
constituted some of the major central nodes in this network. B. In 
papRCC, the most important biological functions of the top deregu- 
lated mlRNAs were: 1) Cancer (p = 9.4E-09-4.15E-02); 2) Endocrine 
System Disorders (p = 9.4E-09-4.54E-02); 3) Reproductive System Dis- 
ease (p = 9.4E-09-5.1 1 E-03); 4) Inflammatory Disease (p = 2.1 7E-06-2.1 7E- 
06) and 5) Inflammatory Response (p = 2.17E-06-1.65E-03). The associ- 
ated functions of the major miRNA network (score = 22) were: 
Endocrine System Disorders, Reproductive System Disease, Cellular 
Development. Tumour protein p53 (TP53) and EPH receptor B6 (EPHB6) 
constituted some of the major central nodes in this network. C. In 
chRCC, the most important biological functions of the top deregulated 
mlRNAs were: 1) Cancer (p = 8.99E-07-4.73E-02); 2) Inflammatory 
Disease (p = 1.64E-06-1.64E-06); 3) Inflammatory Response (p=1.64E- 
06-1.19E-03); 4) Renal and Urological Disease (p = 1 .64E-06-7.96E-03) 
and 5) Reproductive System Disease (p = 3.02E-04-2.39E-02). The 
associated functions of the major miRNA network (score = 22) were: 
Hereditary Disorder, Skeletal and Muscular Disorders, Developmental 
Disorder. Tumour protein p53 (TP53), B-cell CLL/lymphoma 2 (BCL2) 
vascular endothelial growth factor A (VEGFA) and v-myc myelocyto- 
matosis viral oncogene homolog (MYC) constituted some of the major 
central nodes in this network. D. In UT-UC, the most important 
biological functions of the top deregulated mlRNAs were: 1) 
Inflammatory Disease (p = 2.7E-12-3.42E-02); 2) Inflammatory Response 
(p = 2.7E-12-1.75E-05); 3) Renal and Urological Disease (p = 2.7E-12- 
3.35E-09); 4) Cancer (p = 4.67E-10-4.53E-02) and 5) Reproductive System 
Disease (p = 3.82E-06-4.96E-02). The associated functions of the 2 major 
miRNA networks were: 1) Connective Tissue Disorders, Inflammatory 
Disease, Inflammatory Response (score = 25); and 2) Cancer, Reproduc- 
tive System Disease, Renal and Urological Disease (score = 24). Insulin, 
hydrogen peroxide and ribosomal protein SI 5 (RPS15) constituted 
some of the major central nodes in the first network; whereas tumour 
protein p53 (TP53) constituted a central node in the second network. 
doi:1 0.1 371/journal.pone.0091 646.g01 2 

in AB8/ 1 3 cells .significantly reduced the endogenous levels of 
TFCP2L1 protein -21% of total expression; while miR-489 
inhibitors augmented TFCP2L1 levels by reducing the endoge- 
nously expressed miRNA level by -17% (p<0.001; ANOVA) 
(Figure 13). 

Discussion 

RCC subtype classification based on gene expression largely 
reiterates the one based on the histological appearance. However, 
miRNA expression profiles have shown that microRNAs are 
better indicators for distinguishing cancer from normal tissues, and 
can successfully classify even poorly differentiated tumours [30] . In 
the present study we analyzed the global microRNA expression of 
the three most common RCC subtypes along with UT-UC and we 
compared their miRNA profiles with the corresponding profiles 
from normal kidney tissue. The majority of the deregulated 
miRNAs was down-regulated in all RCC subtypes and UT-UC, a 
pattern consistent with miRNA expression profiles in human 
cancers [30]. Our hierarchical clustering also showed that RCC 
and UT-UC tumors could be distinguished perfectly from the 
normal kidney tissue. Given that RCCs and UT-UCs are distinct 
tumour entities, it was not surprising that they formed clearly 
separated clusters. Apart from one chRCC sample which was 
more heterogeneous and shared equal miRNA expression pattern 
with several ccRCC samples, each RCC subtype exhibited a 
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Figure 1 3. Western blotting. A. Normalized luciferase relative light units (RLUs) in AB8/1 3 cell lysates after transfection with sensor constructs. Co- 
transfection of pMIR-REP0RT-TFCP2L1 with either miR-489 mimic or inhibitor resulted in significant reduction or increase in luciferase expression 
versus the negative control, respectively, indicating that miR-489 binds directly on the 3'UTR of TFCP2L1 ("*", p<0.001; "p<0.05, ANOVA). B. Western 
blot of TFCP2L1 from AB8/13 cells after transient transfection with miR-489 miRNA LNA mimics. Inhibitors and the AllStars™ Negative Control 
scrambled sequence LNA. This is a representative of three experiments.The statistical analysis of western blot densitometry results, normalized 
against the Negative Control is also depicted. Values represent the mean ± SEM. Results illustrate the reduction of TFCP2L1 protein levels at the 
presence of miR-489 mimics ( , p = 0.026), while miRNA Inhibitors significantly increased TFCP2L1 levels ( , p = 0.039). 
doi:1 0.1 371 /journal.pone.0091 646.g01 3 



distinct miRNA profile (Figure 3). Furtliermore, in agreement 
with previous reports [5,16,20], ccRCC and papRCC subtypes 
were more closely related when compared with chRCC, reflecting 
their common origin in the proximal tubule of the renal cortex. 

Using a strict filtering approach in the ROC tests, following 
qRT-PCR validation for the top deregulated mlRNAs in each 
tumor entity, we concluded in a group of miRNAs that were also 
validated in an independent patient cohort. Therefore, we suggest 
that these miRNAs could be successfully used as distinguishing 
markers for each RCC subtype or UT-UC. We showed that miR- 
3648, miR-489, miR-638, miR-3656, miR-3687, miR-25-5p, 
miR-21-5p and miR-663b have significant potential diagnostic 
value for ccRCC. Moreover, miR-3648, miR-638, miR-3656, 
miR-3687, miR-25-5p and miR-656 were the best markers for 
papRCC distinction; whereas miR-3687, miR-4284 and miR-141- 
3p can accurately discriminate chRCC over the normal kidney. 
These results were also confirmed using linear regression analysis, 
where those miRNAs manifested both difiFerential expression and 
linear behavior between RCC, UT-UC and the normal kidney. 
However, despite the strict filtering approach that we followed, a 
validation in a larger papRCC, chRCC and UT-UC sample 
number would strengthen the diagnostic power of the latter 
miRNAs. It is particularly important to note that several validated 
miRNAs were suggested as discriminatory markers in more than 
one tumour entities and especially miR-3687 could discriminate 
evenly all RCC subtypes and UT-UC vs. the normal kidney, as 



also shown from our regression analysis. This finding denotes that 
deregulation of miR-3687 seems to be a common manifestation 
during the renal carcinogenic process. Similar miRNA panels were 
recently used by others in order to distinguish various RCC 
subtypes from the normal kidney; however the results are not 
always consistent among them [12,13,16,17,19,20,23,24,31,32]. 
This can be attributed to various factors, such as different 
microarray platforms used, different normalization methods 
followed, difiFerent sample number, etc. Spector et al. just recently 
developed a similar miRNA-based diagnostic assay for the 
classification of ccRCC, papRCC, chRCC and benign oncocyto- 
ma. Following a similar approach, they validated a classifier set of 
24 miRNAs to an independent cohort and they distinguished miR- 
210, miR-195, miR-182, miR-31-3p and miR-551b, among others 
[33]. Interestingly enough, these miRNAs were also included 
among the top 30 deregulated miRNAs in our study (Table S2). 
Apart from the recent study of Izquierdo et al. [34], where the 
authors investigated miRNA expression in progressing and non- 
progressing UT-UC patients, this is the first report to our 
knowledge suggesting a group of miRNAs able to discriminate 
UT-UC from the normal kidney. Despite the small UT-UC 
sample number examined, our data show that miR-3648, miR- 
638, niiR-3656 and miR-3687 were the best molecular markers 
for distinction of UT-UC over the normal kidney (Table 3). We 
also compared the data between UT-UC and RCC subtypes, and 
identified miRNAs that were specific for each tumour entity. 
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Among them 16 miRNAs were specifically down-regulated and 5 
miRNAs were specifically up-regulated in UT-UC (Table S3). 
Approximately 50% of all annotated human miRNA genes are 
located at fi-agUe sites or areas of the genome that are associated 
with cancer which are prone to breakage and rearrangement in 
cancer ceUs [35,36,37,38,39,40,41,42,43,44]. Our data showed a 
partial matching pattern between the reported cytogenetic changes 
and miRNA deregulation in RCC and UT-UC. This is in 
agreement with the recent findings of Chow et al. who 
demonstrated that chromosomal aberrations in ccRCC can 
control the deregulation of miRNAs [12]. 

It was not surprising that IPA revealed that the most significant 
biological functions of the deregulated miRNAs in RCC involved 
cancer and renal and urological disease. Furthermore, the systemic 
inflammatory response has been revealed to contribute to 
tumorigenesis in the last two decades. Key inflammatory 
pathways, such as the NF-kB [45], AP-1 [46] and STATS [47] 
pathways, are implicated in the proliferation, transformation, 
survival, invasion, angiogenesis, metastasis, chemoresistance and 
radioresistance of cancer [48]. Systemic inflammatory response, 
represented by elevated C-reactive protein (CRP), platelet count 
(PC) and erythrocyte sedimentation rate (ESR) was recently shown 
to predict poor survival in patients with RCC [49]. The tumor 
suppressor p53 (TPS 3) and v-myc avian myelocytoniatosis viral 
oncogene homolog (MYC) constituted central nodes in most 
miRNA networks. TP53 plays a central role in tumor prevention. 
As a transcription factor, it mainly exerts its function through 
transcription regulation of its target genes to initiate various 
cellular responses. Recent studies have demonstrated that TP53 
cross-talks with miRNAs at multiple levels; from inducing their 
transcription expression and promoting their maturation, to being 
negatively regulated by them [50]. MYC expression was recently 
shown to be significantiy upregulated in RCC [51] and its 
implication with miRNAs in RCC was also recently reported [52]. 
As a transcription factor, it was recently shown to regulate genes 
involved in cell growth, including those regulating ribosome and 
mitochondrial biogenesis and intermediary metabolism; and 
glycolytic genes, enhancing glycolysis, glutamine uptake and 
metabolism, and lactate production [53]. 

Finally, we isolated the top co-upregulated and top co- 
downregulated mRNA transcripts from 5 independent GEO 
ccRCC datasets, and combining our miRNA results by a 
computational approach, we concluded to the six most anti- 
correlated mRNA/ miRNA parrs the association of which was 
functionally investigated. We showed that the developmental 
transcription factor TFCP2L1 (transcription factor CP2-like 1) 
whose loss was recently suggested to play a role in the loss of 
epithelial differentiation and function in ccRCC, is a true target of 
miR-489, which was among the ccRCC-specific miRNAs in our 
study [54]. 

The current study adds a number of novel deregulated miRNAs 
accompanied by the molecular pathways in which they are 
involved mainly in ccRCC. These miRNAs accurately discrimi- 
nate ccRCC from the normal kidney and can provide the basis for 
a specific targeted therapy. Due to the small sample number of the 
papRCC, chRCC and UT-UC samples that were used in the 
present study, the respective miRNAs should be carefully taken 
into consideration. 

Supporting Information 

Figure SI Pair-wise Hierarchical Glistering (HCl). Pair- 
tissue-wise unsupervised two-way hierarchical clustering with 
Euclidian distance. The log2 fold change in each RCC subtype 



and UT-UC versus the normal kidney tissue was used to construct 
the heat map. miRNA profiling accurately discriminated between 
ccRCC and papRCC (A), ccRCC and chRCC (B), chRCC and 
papRCC (C), ccRCC and UT-UC (D), papRCC and UT-UC (E), 
and chRCC and UT-UC (F). ccRCC, clear ceU renal ceU 
carcinoma; papRCC, papillary renal cell carcinoma; chRCC, 
chromophobe renal cell carcinoma; UT-UC, upper tract urothe- 
lial carcinoma. Red and blue colours show significant up- or 
down-regulation of each miRNA in the tumour versus the normal 
kidney, respectively. 
(TIF) 

Figure S2 Regression analysis. Regression analysis of top 
up-regulated (A) and top down-regulated (B) mlRNAs. miRNAs 
were tested with one-way ANOVA and further linear regression 
was used to fit their expression into a linear model. Left columns 
represent the ANOVA results. Right columns represent the linear 
regression results after sorting of the median values of the 
respective miRNAs in ascending order. 
(DOCX) 

Figure S3 Regression analysis. Regression analysis of top 
up-regulated (A) and top down-regulated (B) mlRNAs. miRNAs 
were tested with one-way ANOVA and further linear regression 
was used to fit their expression into a linear model. Left columns 
represent the ANOVA results. Right columns represent the linear 
regression results after sorting of the median values of the 
respective miRNAs in ascending order. 
(DOCX) 

Figure S4 Regression analysis. Regression analysis of the 
top deregulated miRNAs. One-way ANOVA (left columns) along 
with linear regression (right columns) of the RCC grades with 
respect to miRNA expression was used to fit miRNA expression 
into a linear model. Each regression consists of two curves: one of 
the form y = ax+b (red) and one of the form y = ax3+bx2+cx+d 
(blue). All results presented obtained a p<0.05. The median 
expression values of the mlRNAs were used. Interestingly, miR- 
148b-5p, miR-1910 and miR-200b-5p, miR-210 and miR-3656, 
presented linear behavior in a descending order of the RCC grade; 
whereas miR-3687 and miR-654-5p manifested linear behavior in 
ascending order of the RCC grade. 
(DOC) 

Figure S5 Mviltiple sequence alignments (MSA). Repre- 
sentative MSA and phylogenetic trees of mLR-489 (A) and miR- 
25-5p (B). The majority of the DE miRNAs (78.45%) were highly 
conserved among species. Specifically, 70.68% of the up-regulated 
and 81.22% of the down-regulated miRNAs were phylogenetically 
consei-ved among species. 
(TIF) 

Table SI The table lists the primer sequences that were 
designed to introduce a Spel and a Hindlll restriction 
enzyme sites to be cloned into the pMiR-REPORT 
luciferase vector. 

(DOC) 

Table S2 Top deregulated miRNAs in ccRCC, papRCC, 
chRCC and UT-UC vs. the normal kidney tissue. 

(DOC) 

Table S3 Co-deregulated and tissue specific miRNAs in 
ccRCC, papRCC, chRCC and UT-UC. Co-upregulated and 
tissue specific mlRNAs in RCC subtypes (A). Co-downregulated 
and tissue specific miRNAs in RCC subtypes (B). Co-upregulated 
and tissue specific miRNAs in RCC and UT-UC (C). Co- 
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downregulated and tissue specific miRNAs iii RCC and UT-UC 

(D). 
(DOC) 

Table S4 Reported chromosomal abnormalities associ- 
ated with ccRCC, papRCC, chRCC and UT-UC, as weU 
as some of the well-studied genes located within those 
aberrated regions that are involved in the pathogenesis 
of the tumour subtypes. Many of the significantly deregulated 
miRNAs are located within cytogenetic regions that are signifi- 
cantly aberrated in the difTerent subtypes. In the reported 
chromosomal aberrations the minus sign denotes a loss and the 
plus sign denotes a gain in chromosomal material. Red colour ( | ) 
denotes elevated miRNA expression and blue colour ( J, ) denotes 
reduced miRNA expression. 
(DOC) 
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